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Abstract: We study the spectrum of BPS domain walls within the parameter space 
of N'= 1 U(A^) gauge theories with adjoint matter and a cubic superpotential. Using a 
low energy description obtained by compactifying the theory on x S*^, we examine 
the wall spectrum by combining direct calculations at special points in the parameter 
space with insight drawn from the leading order potential between minimal walls, 
i.e those interpolating between adjacent vacua. We show that the multiplicity of 
composite BPS walls - as characterised by the CFIV index - exhibits discontinuities 
on marginal stability curves within the parameter space of the maximally confining 
branch. The structure of these marginal stability curves for large N appears tied to 
certain singularities within the matrix model description of the confining vacua. 
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1. Introduction 

Four-dimensional J\f~ 1 supersymmetric gauge theories are believed, and in some 
cases have been shown, to exhibit many of the subtle and physically relevant phases 
seen within the standard model. In particular, in the last decade technical advances 
have allowed us to explore properties such as abelian confinement and chiral symme- 
try breaking within specific J\f = 1 theories [1, 2]. However, this progress has generally 
been limited to particular examples, and one may hope that a more global perspec- 
tive on the space of J\f— 1 theories is attainable. In this regard, interesting recent 
work by Ferrari [3] and Cachazo, Seiberg and Witten [4, 5] has provided insight into 
the structure of the quantum parameter space oi Af — 1 theories. This work was 
stimulated by the realisation of Dijkgraaf and Vafa [6] that a matrix model structure 
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apparently underlies the chiral sector of N'= 1 gauge theories, and specifically those 
with adjoint matter. 

For a theory with gauge group U(A^), and a super potential tr>V($) for the 
adjoint chiral superfield the parameter space in question corresponds to the set 
of dimensionless variables, modulo symmetries, that one can construct based on the 
parameters in W($) and the dynamically generated SU(A^) scale Ajv=2- We will 
focus on the simplest nontrivial example, with a cubic superpotential, 

>Vci($) = ^m$2 + i^$3_ 

Classically this theory has two vacua within which the gauge group is unbroken, 
which is the phase on which we will focus. In the limit that g — the infrared 
theory is simply J\f — 1 super Yang- Mills (SYM) and the dependence of >V|„ on 
m/Ajv=2 is fixed by a decoupling relation so that only a single fractional instanton 
can contribute in the confining vacua of the unbroken SU(A^). More generally, for 
finite g one finds a nontrivial sum of fractional instanton contributions, leading to 
an effective superpotential of the form [3] , 

m = I^Aki (l ± [1 - \e''''''^f") , = 0, . . . , AT - 1, (1.2) 

in terms of the dimensionless parameter, 

A = ^^^, (1.3) 

which is a convenient coordinate [3] for the parameter space associated with the 
confining vacua (1.2). 

These vacua lie on the maximally confining branch where the only remaining 
massless field is a decoupled U(l) multiplet. However, this branch connects to oth- 
ers associated with classical vacua for which the gauge group is partially broken. It 
has been conjectured that all these transitions are of two basic types [3]: (1) where 
there are additional massless monopoles, such as the Scibcrg-Witten singularities 
connecting Coulomb and Higgs branches in M= 2 SYM; and (2) where the gluino 
condensate vanishes corresponding to branch points in (1.2). These two cases may 
alternatively be characterised as the result of summation over instanton and respec- 
tively fractional-instanton contributions. In the case of (1.2), the branch points at 
[7, 3] 

Ajfc = (1.4) 

also imply, for N even, the presence of massless monopoles, and thus play a dual 
role. These points were first noted in the latter context, namely as singular points 
connecting Coulomb and confining branches, by Argyres and Douglas for SU(3) [7]. 
Furthermore, in addition to the presence of these singularities connecting branches 
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in different phases, it was shown by Cachazo, Seiberg and Witten [4] that particular 
phases may also exhibit multiple classical limits, connected by smooth transitions 
through strong coupling regions in parameter space. For the cubic model (1.1), these 
smooth transitions always correspond to massless branches with N > 4. 

The presence of a connected set of vacuum branches, in general describing a 
multi-sheeted covering of the quantum parameter space, suggests that it may be 
profitable to study the action of symmetries on this space. To this end, the aim of 
this paper is to go beyond the vacuum structure and explore the spectrum of 1/2- 
BPS states within the parameter space. For the A/'= 1 theories studied here, these 
BPS states are domain walls connecting discrete massive vacua, for which the data 
which can be extracted from the chiral sector of the theory is limited to the tension, 
determined immediately via the vacuum condensates, and also the multiplicity. The 
latter degeneracy is formally determined by the CFIV index [8] in the dimensionally 
reduced 1-l-lD theory obtained by compactifying the worldvolume dimensions of the 
wall on a 2-torus. For our purposes, the CFIV index will be understood to provide 
a definition of the wall multiplicity. 

To determine the multiplicity of BPS states we require a Wilsonian low energy 
description, and since the vacua of interest are confining this requires some defor- 
mation of the original theory. We find that compactifying the theory on a circle of 
radius R, and using the known relation to integrable systems [9, 10], is useful for 
this purpose. This approach also provides a straightforward means of reproducing 
the vacuum structure, as first utilised by Dorey for the M—l* theory [11], and used 
recently by Boels et al. [12] to determine the vacuum condensates for J\f— 1 models 
with adjoint fields, such as the example studied here. 

We explore the multiplicity of certain BPS walls connecting confining vacua as 
a function of the relevant parameter A, and find that the spectrum exhibits discon- 
tinuities on curves of marginal stability (CMS), with certain bound states excluded 
from compact domains in the A-plane. The data describing the multiplicities of BPS 
states thus form nontrivial sections over the parameter space, in a rather close anal- 
ogy to the way BPS particle multiplets form sections, transforming under subgroups 
of SL(2,Z), over the moduli space of J\f—2 gauge theories. While we do not yet see 
evidence for nontrivial quantum symmetries of the latter form, we hope that further 
analysis in this direction will lead to novel constraints. 

The plan of the paper is as follows. In the next section, we describe the compact- 
ification of the theory on R'^ x S*^, which allows a low energy effective superpotential 
to be obtained. This arises from instantons and thus has a natural semiclassical 
interpretation. Using this effective theory, we reproduce the condensates (1.2) in the 
maximally confining vacua. In sections 3, we turn to BPS walls and study the struc- 
ture of the central charges as a function of A, and exhibit a class of marginal stability 
curves for composite states. The structure of these curve simplifies for large N and 
has an interesting connection with certain singularities within the relevant Dijkgraaf- 
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Vafa matrix model, which we also describe. In Section 4, we turn to the dynamical 
question of whether these marginal stability curves do indeed signify discontinuities 
in the spectrum, and verify that this is the case by first deducing the spectrum at 
a special point - A = - and then constructing the leading order inter-wall poten- 
tial as a function of A. This leads to the conclusion that specific composite states 
are removed from the spectrum within a compact domain of the A-plane. Section 6 
contains some concluding remarks, along with further comments on the embedding 
of this theory within M^l* SYM. 

2. Compact ificat ion and the Vacuum Structure 

The class of theories we will focus on here contains an M— 1 vector multiplet with 
gauge group U(iV) and an adjoint chiral multiplet with scalar component 0. The 
cubic superpotential 

tr>Vci($) = ^mtr$^ + \gir^^ (2.1) 

leads to a set of classical vacua, distinguished by the distribution of the eigenvalues 
of between the two minima, 

TTX 

= 0, and (/?(2) = - — . (2.2) 

With A^i eigenvalues equal to ip{i) and N-2 equal to ^{2), satisfying A^i + = iV, 
the gauge group is broken to U(A^i) xU(A'2). We will be concerned primarily with 
those vacua where the gauge group is classically unbroken, i.e. with {Ni,N2) equal 
to (iV, 0) and (0, N). Note that the overall U(l) factor of the gauge group is central 
and thus decouples as all fields are correspondingly uncharged. 

Our aim in this section will be to construct a low energy effective superpotential 
suitable for use in extracting the BPS spectrum. This will necessarily involve a 
deformation of the theory. However, to first deduce the quantum vacuum structure 
on the confining branch, we can proceed more directly. Treating this system as a 
perturbation of Af= 2 SYM, the vacuum structure follows from the Seiberg-Witten 
solution, i.e. given Wd, we can write 

Weff = (trWd(pI + ^)) = NWM + ^}vS\p){tT^^), (2.3) 

where $ transforms in the adjoint of SU(iV), and so we have made use of the relations 
(tr$^"'"'"^) = 0, while p = (tr($))/A^. For confining vacua, the maximal degeneration 
of the SU(A^) Seiberg-Witten curve [13, 14] leads to the non-vanishing condensate [15] 
(tr<l>^) = 2A^cjfcA^^2 with 00^ = exp(27riA;/iV) an A^*^-root of unity. On substitution 
into (2.3), one finds [3, 16] 

Wk{p) = N [WM + m{p)u;kAl=2\ , A; = 0, . . . , TV - 1, (2.4) 
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where m(p) = W"i(p) is the effective mass term for all the scalar modes, including 
the trace component p, at weak gauge coupling. For the cubic case this result, first 
obtained in [3, 16] , has the simple interpretation of contributing the gluino condensate 
from the confining SU(A'") factor. More generally, one expects a nontrivial (but finite) 
sum of such fractional instanton contributions* . Integrating out the trace component 
p, we recover the confining vacua (1.2) noted earlier. 

When the trace component p is light relative to the other modes - and indeed it 
becomes massless when A = - the effective superpotential (2.4) may provide a reli- 
able low energy description. However, since (2.4) depends explicitly on the vacimm, 
we cannot use it to study wall configurations which interpolate between different 
vacua. These configurations must necessarily excite the gauge modes responsible for 
gluino condensation, and we will need to include them in any tractable low energy 
description. Moreover, since the vacua of interest are confining, some deformation 
of the theory is required. To this end, if we wish to avoid adding additional fields, 
a natural approach is to compactify the theory on a circle of radius R. A Wilson 
line for the gauge field then provides a massless scalar, which may be rotated to 
the Cartan subalgebra Pa = J^i Aa, a = 1 . . . N, and this will generically Higgs the 
gauge group to its maximal torus. The resulting photons can be dualised to periodic 
scalars cr" [18], and supersymmetry then ensures that the complex combination, 

qa^ Pa + raa, a^l,...,N, (2.5) 

where r is the complex gauge coupling, is the lowest component of a (classically 
massless) chiral supcrficld which characterises the gauge sector of the dimensionally 
reduced low energy theory. The additional adjoint scalars reduce straightforwardly, 
and the effective 3D theory is a Wess-Zumino model with two adjoint chiral fields and 
a classical superpotential given by (2.1). We now turn to possible nonperturbative 
corrections. 

2.1 Semiclassical interpretation of Wes 

When the compactification radius R is sufficiently small, AR <€. 1, the system is 
weakly coupled and thus any nonperturbative corrections to W should have a semi- 
classical interpretation in terms of instantons. Therefore, we can parametrise the 
superpotential in the form 

Weff = >Vel(p„) + mnstiqa,Pa), (2.6) 

where Wciass(Po) is the classical superpotential (2.1) and we have now chosen to 
denote the eigenvalues of the adjoint field $ as {pa}, a — 1, . . . , N. 

*In certain cases, there may also be additional instanton contributions, whose form determines 
the ultraviolet completion of the theory. See [17] for a nice discussion of this issue. 
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The form of the instanton generated superpotential in these theories was first 
explored by Polyakov [18], and for theories with M=2 SUSY by Affleck, Harvey and 
Witten [19]. For pure M—1 SYM, there are nontrivial 1-instanton configurations, 
corresponding to A" — 1 BPS monopoles ahgned along simple roots and, for finite R, 
an additional 'Kaluza-Klein'-monopole wrapped around the compact direction [20], 
which contribute to the superpotential [21, 22, 23, 24, 25, 11, 26, 27, 28]. In the 
presence of additional adjoint chiral fields, these configurations would normally have 
too many fermionic zero modes to contribute to W, but if these fields are classically 
massive, the additional zero modes can be absorbed by the mass terms. Accounting 
for the monopole-antimonopole background [18] which resums the instanton vertices 
to a field-dependent superpotential, we can write down its generic form as follows 
[26, 27]t, 

1 ^ 

>Vi„,t = 2^k2 ^[^{Pa) + m(p„+i)]e^'^-«'^+\ m{pa) = m + 2^p„, (2.7) 

a=l 

in which m{pa) is the classical mass term for ipa given by expanding Wd about a given 
vacuum (pa). Note that Qa and Qa in (2.5) are related by a constant shift depending 
on the coupling r. The structure of the mass insertions follows from the choice of 
basis for the scalars which corresponds to taking projections along the fundamental 
weights. The monopoles in contrast are aligned along the simple roots, and thus are 
embedded along two fundamental weight directions. 

In the expression (2.7), we have slightly extended the result of Davies et al. 
[26, 27], having made the replacement (pa) Pa = (Pa) + ^Pa, so that a genuine 
interaction term is introduced. We have not derived the appropriate instanton vertex 
explicitly, as it can be justified in a different manner to be discussed shortly. In 
particular, for later use, it is convenient to recall another approach to this problem 
which starts from the A/'=2 hmit and uses the relation to integrable systems [11]. 

2.2 The effective Toda superpotential 

At a formal level, the calculation of Wcfr is dramatically simplified by recalling the 
relation between the Seiberg- Witten solution for JV=2 SYM and certain (complexi- 
fied) integrable systems [9, 10]. More precisely, written in terms of the natural gauge 
invariant coordinates on the N'—2 moduli space, 

ii„=itr(/)", (2.8) 

the superpotential takes the form 

trW ^mu2 + gus, (2.9) 

'^For simplicity, we will write the superpotential using a 3-l-lD normalisation in what follows, so 
that WsD = 27ri?Weff. 
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where we have now dropped the subscript as one can argue that this form is exact 
within the compactified theory on x 5^. To do this, one treats this superpotential 
as a perturbation of the (reduced) jV= 4 theory, and interprets the coordinates 
Un — Un{x) as coustraiued variables depending on the underlying Seiberg-Witten 
curve. To make use of this expression one must find the corresponding unconstrained 
variables parametrising the Jacobian of the curve. 

This problem is elegantly solved via the correspondence with integrable systems, 
and in this case the affine Toda lattice [9] . This relation is perhaps most transparent 
on noting that the Seiberg-Witten curve for U{N) can be written in the form, 

A2JV ^ 

Pn{x, Un) = z + Pm{x, Ur,) =x'' + ^-^"^"^ (2.10) 



z 



n=l 



where s„ and Un are related by the recursion relation, rsr + Y^i=o ^r-iUi — 0, with 
Uq — O and Sq = 1- This curve is identifiable as the spectral curve for the affine Toda 
lattice, namely 

det(xI-L(z)) = 0, (2.11) 
where the Lax matrix is given by 



Pi 
1 




Ajv=2 e«i-«2 



P2 

1 



A 



N=2e 

P3 








(2.12) 



in terms of the canonical variables {pa, Qa}- Thus the unconstrained variables on the 
Jacobian of the curve are precisely the (complex) coordinates and momenta {qa,Pa} 
for the affine Toda lattice. The conserved quantities, to be identified with 
given by 

(2.13) 



Un = -trL", 
n 



which allows us to identify the momenta Pa with the eigenvalues of the adjoint scalar 
Pa — (fa), explaining the notation introduced above. 

With this parametrisation, the superpotential (2.9) takes the form 



N 



N=2^ 



a=l 



N 



Y,Pl + ^^l=2{Pa+Pa+l)e'^ 



-la+l 



o=l 



(2.14) 

which, on rewriting it in the form (2.6), agrees with the expected 1-instanton cor- 
rected potential*. 

*Note that if we were to consider a classical superpotential of degree four or higher, e.g. tr4> , 
then the semiclassical interpretation would require the contribution of certain multi-instanton cor- 
rections. It would be interesting to see precisely how the additional zero modes are lifted in such 
cases. 
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This latter approach for calculating Wcs, making use of the connection to inte- 
grable systems, was first used by Dorey for the J\f = 1* theory, and has since been 
applied to other theories including pure Af— 1 SYM [28, 27] and more general defor- 
mations Wci [29, 30]. A wider class of U(A^) examples involving polynomial superpo- 
tentials was also recently studied in this way by Boels et al. [12], and subsequently 
extended to other gauge groups [31] . 

2.3 Quantum vacuum structure 

Extremising the superpotential (2.14), we find the following equations, 

PaiVPa + 1) = = -viXa + Xa+l) (2.15) 
{I + V{Pa + Pa+l)) Xa = + V{Pa-l + Pa)) Xa-U (2-16) 

where 77 = g/m and Xa — expiQa — Qa+i), and we have imposed the decoupling 
condition Qn+i = <ii for the overall U(l) factor. Since we are interested in vacua with 
a classically unbroken gauge group, we solve (2.16) by setting 

?a = 5, Pa=P, Va=l,...,A^. (2.17) 
The resulting 2A'" vacua are given by 

«i = = ^Pt = Y{^^ Vl - Axfe) , (2.18) 

a 

with 

= e^^^^/^, = 0, . . . , iV - 1, (2.19) 

where -Ui, and thus p, is single valued on a double-sheeted cover of the A-plane, and 
the two sheets arc distinguished by the corresponding classical vacuum expectation 
value for (f). One may note here that the intuitive correspondence between massive 
vacua and (complex) mechanical equilibria of the integrable system, that is so ap- 
parent for SU(iV) [11], appears to have been lost here due to the nontrivial vacuum 
expectation value for p, so that the 'equilibria' now have nonzero coordinate momen- 
tum. However, one finds that these vacua still correspond to stationary points where 
the angular momentum within the Jacobian vanishes [32, 33]. In other words, there 
is a canonical transformation within the integrable system to an action-angle basis 
for which the massive vacua still map to mechanical equilibria. 
The 2A'' extrema of the superpotential are then given by, 

27V 

>V± = ^A^,^, (1 T [1 - \xuf'^) , = 0, . . . , TV - 1, (2.20) 

and these results are fully consistent with those recently obtained using Seiberg- 
Witten and matrix model techniques [3, 16] as discussed above. 
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Figure 1: Confining vacuum structure in the W-plane for TV = 14 and A <C 1, exhibiting the 
action of the cyclic 'Z2jv symmetry', and illustrating the interpolating profile of a /c-wall. 



For A ^ 1 the vacua fall on two approximately circular curves in the W-plane, 
as exhibited in Fig. 1. For the vacua corresponding to the classical vacuum at = 0, 
the asymptotics near A = have the form 



>V+(A ^ 0) ^ TVA^^^e^-^/^ i 1 + O i ( ^ , (2.21) 




where the leading term corresponds to the pure N'= 1 SYM limit, and the subleading 
terms, proportional to powers of e^'^*'^("*)/^ , correspond to higher order fractional 
instantons. 

The plot also exhibits a discrete cyclic symmetry which rotates the vacua via its 
action on A, 

A ^ e'^^/^A, (2.22) 

which we interpret as induced by SL(2,Z) translations on the bare gauge coupling 
T — > T + 1. The corresponding action on the vacua within each branch is then 

- ^."=+1, (2-23) 

while the rotation of A around one of the branch points induces the action [3] 

W,^ ^ (2.24) 

With a suitable choice of branch, the action (2.22) rotates Wfc through all 2N vacua 
on both branches, and we can interpret this 'Z2iv cyclic symmetry' as a nontrivial 
extension of the Z^r nonanomalous discrete subgroup of the classical U(1)r symme- 
try, which rotates the vacua of pure SYM. The extension arises directly from the 
double-sheeted structure of the confining branch over the quantum parameter space, 
coordinatised by A. 



Note that this structure does not, however, extend globally over the parameter 
space. In particular, the asymptotics near X = oo are 

WiiX - oo) ~ ±^-^NgA%^,e'-^/- (' ^ ^ (a^)')) ' ^''"'^ 

and A is no longer a good coordinate on the parameter space - the leading dependence 
is instead on a/A. 

We have focused on the maximally confining vacua, but for completeness we note 
that for generic N there are in addition a large number of branches with massless 
vacua which map to classical vacuum solutions with a broken gauge group in the limit 
that AAr=2 0. The cubic superpotential only exhibits classical vacua where the 
the eigenvalues {pa} are split into two groups, and thus there is a single (nontrivial) 
massless U(l) factor. The structure of these vacua was elaborated for A?" < 6 in [4], 
and shown to include branches where smooth transitions in parameter space between 
distinct classical limits were possible. These vacua decouple in the hmit that A — > 0, 
where the theory reduces in the infrared to pure SYM. However, in the opposite 
limit, A oo, certain massless vacua survive and in the U(3) case, for example, 
apparently asymptote to the Argyres-Douglas conformal points [7] . 

3. BPS Wall Kinematics 

With a clear picture of the confining vacuum structure in the W-plane, we can now 
turn to the question of the spectrum of BPS states, namely domain walls interpolating 
between distinct vacua. Given that the vacua described above are distinct, such 
configurations are necessarily present on topological grounds. However, the number 
of corresponding BPS states is a dynamical question that wc will come to shortly. 
We first outline the kinematic structure imposed by the J\f= 1 superalgebra. 

3.1 Central Charges 

The vacua labelled by Wk allow us to isolate two sets of walls, characterised by 
whether or not both vacua at spatial plus and minus infinity are on the same branch. 
If they are, then such walls are present only at the quantum level. However, if the 
asymptotic vacua are on different branches then the corresponding walls are visible 
classically. More precisely, as noted in [3], this characterisation is useful for A <^ 1 
while, as we will see, by varying A one observes various discontinuities in the spectrum 
of BPS walls, with the underlying reason being the presence of curves of marginal 
stability in the parameter space coordinatised by A. 

We begin by discussing the central charge structure, and to this end it is useful 
to introduce the following (condensed) notation for the vacua, 

27V 

m = ^Aki (1 - [1 - Axfe]'/') , A; = 0, . . . , 2Ar - 1, (3.1) 
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which incorporates both branches, with the imphcit understanding that the vacua 
for k = N, . . . , 2N — 1 he on the second branch. This notation is consistent with the 
discrete symmetry noted above. It is then convenient to introduce the (SU(A^) 
root-valued) topological charges, 

^fe) = S'- (3-2) 

in terms of which the central charges for walls interpolating between vacua labelled 
by i and j phase-units respectively are given by [34] 



(3.3) 



The second term here is an anomaly, which vanishes in the present case due to the 
J\f—2 matter content. Consequently, we have 

Zij = 24^.)Wfe = 2 {Wj - Wi) , (3.4) 

or more explicitly, 

while the corresponding wall tensions are 

Tjk — \Zjk\- (3.6) 

Note that T^k+N goes to zero at the branch points in the A- plane [3], corresponding 
to the fact that the vacua labelled by Wfc and Wk+N collide at these points. 

3.2 Marginal Stability 

We now specialise to walls interpolating between two vacua on the same branch, 
and for definiteness restrict to walls associated with the central charges Z_kk- For 
sufficiently small A, such configurations are naturally interpreted as bound states of 
2k minimal 1-walls with charges which interpolate between adjacent vacua. 

For this reason, we will generally focus on the simplest bound state Z_ii composed 
of two 1-walls. Supersymmctry demands that such putative BPS bound states, when 
present, are at least marginally bound. The corresponding submanifolds on which 

r_n = T_io + roi (3.7) 

are of co-dimension one in the parameter space, and allow for possible discontinuities 
in the spectrum, where BPS bound states may delocalise and leave the spectrum. 
The position of these curves of marginal stabihty (CMS) in parameter space is fixed 
by supersymmetry and therefore is purely kinematic. It is this question that we will 
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Figure 2: A plot of the CMS curves in the A-plane for the 2-wall composite, and several values 
of N. For N ^ oo the limiting form of the curve is a circle described by (3.14) 

address in the remainder of this section, while the subsequent dynamical issue of 
whether discontinuities in the spectrum do actually occur on these submanifolds will 

be addressed subsequently. 

A convenient characterisation of the submanifolds on which (3.7) holds is that 
the relative phase of the two constituent central charges vanishes, 

c^(A)|cMS = 0, (3.8) 

where 

^ (3.9) 

1^-10^01 1 

Using the explicit form for the central charges (3.5), the condition (3.8) imphes the 
following two (real) constraints: 

Im(A(A)A(A)) = 0, sgn(A(A)A(A)) < 0, (3.10) 

where 

fkiX) ^ 1 - (l - Si^e-'^/^sin . (3.11) 

The CMS defined by the constraints in (3.10) is exhibited in Fig. 2 for several values 
of A^; it defines a closed curve for N > 6, and has a simple limiting form for large N 
that we will now describe in more detail. 

As follows from (3.5), the expression (3.9) for the relative phase contains branch 
points at A = {e~^'^'/^, 1, e^'^*/^} internal to the A-plane, and this branch structure 
is inherited by the equation (3.10) defining the CMS. Although having physical sig- 
nificance in signalling the presence of singular points connecting this branch with 
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others associated with massless vacua, these points are not crucial here in that the 
spectrum of finite tension walls is at most discontinuous on curves of co-dimension 
one, and so the co-dimension two branch points can always be avoided. 

To make use of this fact, we can simplify the structure of the CMS by taking N 
large, in a scaling limit such that 



whereupon at leading nontrivial order the branch points lie along a line, ~ 1 ± 
2mk/N. In this regime it is straightforward to evaluate the relative phase oj{\) for 
the 2^-11 wall bound state. We obtain the following simple expression. 



which reduces to a; = 27r/A^ in the limit A — > of pure SYM, and more generally 
is a harmonic function of A, with a corresponding pole at the branch point Aq = 1 
associated with the intermediate vacuum. 

Prom this expression, we find that at large N the phase uj vanishes on a circle in 
the complex A-plane described by 



Note that since the curve only touches the line A = 1 + ia at the point a = 0, it 
is consistent to track the curve all the way to A = 1 since it avoids the additional 
branch points at A_i and Ai. There are apparently two distinguished points on the 
CMS, at A = 1 and A = 2/3, where the curve intersects the real axis. The former 
corresponds to the branch point for the intermediate vacuum, while the latter has an 
interesting interpretation within the Dijkgraaf-Vafa matrix model [6], that we will 
explore in more detail below. 

3.3 The matrix model 

The particularly simple large- structure for the CMS determined above hints at a 
more transparent interpretation. One point of view that we will provide some evi- 
dence for below is that the onset of marginal stability in this regime is characterised 
by a qualitative change of the intermediate vacuum. For the composite state consid- 
ered above, the relevant vacuum has a critical point at A = 1, and presumably also 
intersects with massless vacuum branches at other points whose precise locations 
depend on N [3, 16]. 

The fact that the CMS passes through the critical point at A = 1 is then not 
too surprising as the intermediate vacuum at this point will have additional massless 
excitations. The second intersection of the CMS with the real axis at A = 2/3 is 
less immediately attributable to any pathology of the vacuum. However, it is a 



Nuj[\) = constant. 



(3.12) 




(3.13) 



CMS: |6A-5p = l. 



(3.14) 



distinguished point within the geometric description of the vacua in the jV=2 hmit, 
and thus also within the matrix model picture. 

To illustrate this point, we first recall some of the relevant aspects of the Dijkgraaf- 
Vafa matrix model. In particular, for the theory at hand, there now exist some rather 
general arguments [6, 35, 16] implying that specific aspects of the chiral sector of the 
theory are captured by a holomorphic matrix model of the form 

^-M^Hsys^ = J [dX]e-f ^--^^i^^), (3.15) 

where X is an M x M matrix, and we take M — > oo to isolate the planar sector. 
S serves as a loop counting parameter in the planar limit of the matrix model, but 
ultimately, one makes the identification S — —-^^trW^Wa within the gauge theory. 
The large-M saddle point is characterised by the condensation of the eigenvalues of 
X into a set of cuts [ak,bk], each corresponding classically to one of the roots of 
W^i(a;) = 0. The saddle point condition is conveniently expressed in terms of the 
force on a test eigenvalue, 

y^{x) =W',,-2Snm where TZ^ = ^ {Tr . (3.16) 

M \ X — X I 

TZja is the trace of the resolvent whose discontinuity across the cuts determines the 
eigenvalue distribution (see e.g. [36]). We are concerned here with 1-cut solutions, 
corresponding to classically unbroken U(A'^) vacua, and in terms of i/m the (large M) 
saddle point requires [6, 3], 

yl = {yV'Mf - /i(^) = 9'{x - x.Yix -a){x- b) (3.17) 

where fi is a polynomial of degree one, and a and b are the cndpoints of the cut, 
with a double point corresponding to the degenerate second cut. This constraint 
describes a degenerate hyper-elliptic curve, and is equivalent to the appropriate de- 
generation of the U(A'") Seiberg-Witten curve [13, 14]. 

For this saddle point configuration, the gauge theory superpotential is deter- 
mined by matrix model prepotential as [6] 

]A^^ = NUb - {t + k)UA, (3.18) 

in the 1-cut sector, where the periods (11^,11^) are given in terms of the resolvent 
TZra via 

= 2mS = jf y^, Ub = 2m^ = ^ y^, (3.19) 

using a symplectic basis of (compact) A and (noncompact) B cycles for the surface 
(3.17), as shown in Fig. 3. 

Making use of (3.17), the integrals in (3.19) may be expressed in terms of ele- 
mentary functions. Evaluating IIb in the limit Aq — > oo, and dropping an irrelevant 
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constant, one finds that the contribution from 11^ in (3.18) serves to reconstruct the 
dependence on the dynamical scale Ajv=2- With p — p{S) denoting the mid-point of 
the cut, one finds that [6, 3, 16, 37], 



iV=2 



N 



(3.20) 



This result - a generahsation of the Veneziano-Yankielowicz [38] super potential - is 
identical^ to the answer one obtains from (2.4) by integrating in S via a Legendre 
transform with respect to InA^g [16], and then integrating out p in favour of S. 
On extremising VVui{S), and with the implicit dependence on the scale A^=2, the 
on-shell superpotential takes the form. 



where 



dU 



B 



k 



(3.21) 



(3.22) 



OUa N 

is the modular parameter of the noncompact elliptic curve, indicating that it is an 
A'^-fold cover of the bare curve. Of course, a renormalisation group-invariant version 
of this statement is that Wm|it depends on N only via the (complexified) dynamical 
scale 



A^=2 exp 



27riA;^ 



Apexp 



27ri(r + k) 



(3.23) 



N J ^ - y TV 

This explanation of the remarkably simple A^-dependence of the vacuum values of the 
superpotential, and consequently of the central charges is one of the virtues of the 
matrix model approach, since the matrix model itself is A'"-independent, and results 
from the planar saturation of the chiral sector [6]. In particular, using the topological 

^The on-shell equivalence between the 'equilibria' of the integrable Toda system and the matrix 
model saddle points has recently been verified more generally [32, 33]. 
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h=2/3 



Figure 4: A schematic illustration of the degeneration of two specific cycles within the generic 
J\f= 2 U(iV) curve. We can interpret the point A = 1 as corresponding to the degeneration of one 
of the cycles of the generic curve, while the point A = 2/3 corresponds to the degeneration of two 
such cycles. In contrast, the 1-cut saddle-point for the matrix model is always characterised by one 
non-degenerate cycle and a double point at cc*. As shown, at A = 1 the double point sits at the 
midpoint of the remaining cut, while for A = 2/3 it sits at one end. 

charges n^j^-^ = Sj — 5^, the central charges can be compactly written in terms of the 
S-period as follows, 



Geometrically, the 1-cut configuration is unique, with the nontrivial vacuum struc- 
ture arising via extremisation of (3.18) for a given value of k. However, for the 
purpose of interpreting domain walls in this picture, it is useful to note that the 
rotation of the cutoff Aq — > AqB^'^* induces the monodromy 11^ — > 11^ — 11^ [37], 
and thus if we choose to fix the phase of the regulator Aq, we can view the process 
of passing through a wall as corresponding to a rotation of the cut as exhibited in 
Fig. 3, with the corresponding rearrangement of the eigenvalues^. 

With this formalism in hand, we now reconsider the near-CMS regime, and in 
particular the intercept points A = 1 and A = 2/3. As follows from (3.17), both these 
points correspond to degenerations of the 2-cut solution as shown in Fig. 4. In the 
plot, for illustration, we show some of the additional cycles relevant within a 2-cut 
degeneration of the Seiberg-Witten curve. However, in the 1-cut solution we consider 
here, the only remnant of these degenerations is that at A = 1, the double point at 
X = lies at the mid-point of the cut, while at A = 2/3 the double point lies at one 
end. The degenerating cycles do not play a direct role in the matrix model solution. 

The point A = 2/3 is, however, known to correspond to a critical point within 
the 1/M-cxpansion of the matrix model [3, 36], i.e. of the prepotential. This suggests 

^In the context of 'classical walls', where the cuts corresponding to the two vacua have distinct 
classical limits, it was noted by Dijkgraaf and Vafa that walls are apparently dual to eigenvalue 
tunnelling [6]. Such an 'instanton-like' interpretation seems more problematic here as the relevant 
vacua are classically degenerate. 




(3.24) 
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that it should be possible to see consequences of the degeneration within (3.20). At 
A = 1, this is of course reflected in the branch point for VVo(A). Evidence of the 
degeneration at A = 2/3 is less apparent, but can be uncovered if we define the 
following '/3-function' for the period of the curve (3.17), 



dS 



(3.25) 



where the period (3.22) is given (off-shell) by t{S) — dgT in terms of the prepotential. 
Evaluating /3(A) in the relevant A; = intermediate vacuum, we find 

= 27vl^, (3.26) 

which vanishes for A = 1, and diverges for A = 2/3. This direct connection between 
degenerations of the geometry and the position of the CMS for large is rather 
suggestive, but without a clearer picture of the physical excitations associated with 
>Vefr(5') it is difficult to fully interpret the consequences of the latter singularity", 
and we will not pursue this further here. 

In the next section, we turn to the dynamical question of deducing the multi- 
plicity of states, and indeed whether discontinuities do indeed arise on crossing the 
CMS curves discussed above. 



4. BPS Wall Spectrum 

The multiplicity of 1/2-BPS multiplets mj\f=2 theories in 1+lD with fixed boundary 
conditions is given by the CFIV index, i/jk, formally defined as [8] 

ujk = Tr,kF{-lf, (4.1) 

where F is the fermion number of the corresponding state. The crucial property 
of Ujk, following from its definition, is its stability under various deformations of 
the theory. In particular, it was shown in [8] that it is stable under (nonsingular) 
variations of the D-terms, such as the Kahler potential, which are unconstrained 
by supersymmetry. This makes calculation of i^jk tractable as we can deform the 
theory, for example via compactification as here, and the result can be shown to be 
independent of the compactification radius R which enters the Kahler potential. 

Moreover, ujk is also stable under certain variations of the F-terms. For a Wess- 
Zumino model, characterised by a Kahler potential /C(X", X") and a superpotential 
W(X"), as we have here this can be understood via a geometric reformulation of Uj^ 

"As noted recently [39], another quantity with the same characteristics is the formal 'mass term' 
VV"g{S) which, however, is again difficult to interpret without knowledge of the Kahler potential. 
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[40]. Specifically, an important consequence of the Bogomol'nyi equation for 1/2- 
BPS walls is that the superpotential must trace out a straight line in the W-plane 
[41, 42]; i.e. 

>V(X") = (1 - t)Wo + tWi, (4.2) 

where t G [0, 1] is a fiducial variable parametrising the transverse coordinate to the 
wall, and Wo and Wi are the corresponding vacua between which the wall interpo- 
lates. The set of solutions to the hnearised Bogomol'nyi equation about each vacuum 
defines a cycle in field space, and one can show [40] that by comparing the two 
cycles Aj and at a given value oi t & [0, 1], the number of solutions is given 
precisely by the intersection number, 

ujk = A, o Afe. (4.3) 

The topological character of the intersection number then implies that i/jk is stable 
under deformations of W which do not lead to additional vacua crossing the straight 
line trajectory between Wo and Wi, as it is only in this case that the intersection 
number varies according to a Picard-Lefschetz monodromy [40]. 

Thus, to prove the existence of 1/2-BPS states it is sufficient to verify the pres- 
ence of smooth profiles in field space for which the interpolating trajectory in the 
punctured W-plane - with additional vacua excised - is homotopic to the straight 
line (4.2). With this argument in mind, and given our knowledge of the structure of 
the CMS curves, our strategy will be first to determine the spectrum at a convenient 
point - we will consider the pure SYM regime A — > - and then to study whether 
this multiplicity is preserved on crossing the CMS. For the latter test, we construct 
the leading order interaction potential, as a function of A, between the constituents 
of the 2-wall composite discussed in the previous section. 

4.1 Multiplicity for A 1 

We will first consider the spectrum at a special point A = 0, where the infrared 
theory reduces to pure A/'= 1 SYM. BPS walls in this theory have been studied from 
several points of view [43, 44, 45, 46, 47, 48]. Within the low energy description 
on X 5"^, the effective superpotential reduces in the A — > limit to the following 
(complexified) affine-Toda form [21, 22, 23, 24, 25, 11, 26, 27, 28], 

N 

Weff(A = 0) = A^=i (4.4) 

a=l 

subject to the constraint that Yl^a ~ 1- Using the cyclic symmetry, we choose the 
initial vacuum to he Xa — 1 and the second vacuum to be Xa — e^'^*'^/^, with k < N. 
Within the wall, the winding number of each field must then vary from zero to k/N 
mod Z, and moreover the constraint f|xa = 1 implies that at any particular point 
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along the trajectory (4.2) the winding numbers of the fields Xa must sum to zero. 
Thus, at any non- vacuum point, they cannot all be equal. To minimise gradient 
energy, there are only two allowed winding numbers: k/N and k/N — 1 [49], and 
this leads us to the following resolution of the constraint in terms of a new variable 

y [50], 

Xa=l,...,k = Xa=k+1,...,N = V''' ■ (4.5) 

The relation (4.2) then imphes that for the corresponding fc-walls 

ky^-k + (AT - k)y-'' = N{1 -t) + Nte^'"'^'^. (4.6) 

For N — 2 (and thus k — 1), this relation can be straightforwardly inverted to obtain 
the interpolating trajectories. However, for generic N, we can instead simply verify 
the existence of BPS solutions by following an approach used in the context of the 
J\f— 2 affine-Toda model in 1 + lD** [50]. Note that since y is a complex variable, 
there are two linearised solutions about each vacuum, but at most one combination 
of these can link to form a BPS trajectory. To verify that such a trajectory exists 
we follow the argument described above [40, 50]. In particular, with the ansatz 
y(t) = e2^^*/^, we find 

\Wes(X ^ 0)\ ^ \N - k + ke^'^'^l < N, (4.7) 

where the final equality holds only at the initial and final vacua. Therefore, since 
the remaining vacua lie on a circle of radius N, this suffices to prove the existence 
of the corresponding BPS multiplet [50]. Moreover, accounting for the permutation 
symmetry in the identification of y, the full multiplicity is given by 

m{x^o)^(^^y (4.8) 

This result, conveniently interpreted as a multiplet transforming in the k^^ fundamen- 
tal representation of SU(iV), is consistent with other calculations of the degeneracy 
of BPS walls in pure J\f^ 1 SYM [53, 49]. 

4.2 Inter-wall potential for generic A 

To explore how the multiplicity of BPS states varies as we move in parameter space, it 
is useful to follow a somewhat different approach. It should be clear from the analysis 

**As an aside, bosonic affine-Toda theory in I+ID with an imaginary couphng is also known to 
exhibit sohtons [51]. In general, these solutions satisfy a second-order equation with the potential 
(4.4), and so are not related to the configurations we are concerned with here. However, the case 

= 2 is an exception where the sine-Gordon soliton is (up to certain rescalings) a BPS solution if 
we use the classical Kahler metric. One may understand this via noting that affine-Toda solitons 
also satisfy a set of first-order (Backlund) equations [52] which reduce precisely to the Bogomol'nyi 
equations for A'' = 2. 
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above that the spectrum of 1-walls will be stable for |A| < 1 where the structure of 
neighbouring vacua does not change qualitatively. However, it is less apparent that 
the multiplicity of higher A;-walls is stable in this region of the parameter space. 
We can determine possible discontinuities by treating /c-walls as bound states of k 
1-walls and calculating the inter- wall potential as a function of A. In practice, it will 
be sufficient to focus on the simplest case of a 2-wall bound state, corresponding to 
the example considered in the previous section. 

To proceed, we note that on general grounds the only global symmetries broken 
by the wall configurations arc super-translations. We then infer that the only bosonic 
moduli of a system of two asymptotically separated 1-walls will be their respective 
centre-of-mass positions. It is an implicit assumption that any other worldvolume 
fields are massive and can be ignored for the purpose of considering the ground state. 
The bosonic moduh space of the constituent system is then R^, and we can always 
decouple the overall translational mode, and thus the relative moduli space is one- 
dimensional, parametrised by a real field r. The problem at hand then reduces to 
computing the potential induced on this space when the constituent 1-walls are at a 
finite separation. 

For the region |A| < 1, it convenient to further restrict the problem by considering 
the regime where RAn=2 ^ 1 which means that the adjoint scalar fields Pa are 
(generically) of higher mass than the gauge modes Qa and can be integrated out, 
leading to the reduced superpotential. 




which is now a function only of Xa- However, the presence of branch points implies 
that this reduced system is valid only for sufficiently small A, for which the field 
profiles remain well away from the cuts. 

As noted above, we focus on 2-walls and fix the two vacua to be = ^-2mjN 
and Section 3. The putative 2-walls in this sector, counted by 

can be viewed as bound states formed from the 1-walls for which z/_io = and 
t'oi = A^. A schematic illustration of how the phases of the Xa fields must vary over 
the wall profile is shown in Fig. 5. 

The leading order potential between the constituents can be determined by ex- 
panding the tension of wall two in the background of wall one (see also [54, 55]). 
More precisely, we can write the tension of the second wall T2 — 2|W(i) — W(o)| as 

T2(r) = T2 - 2Re [e-^5^io)(;2o)5(oi)(^o + r)a„a;,>V(o)] + • • • , (4.10) 

where the second term is the leading order interaction potential V{t) between the 
constituent 1-walls positioned at z — Zq and z — Zo + r respectively (see Fig. 5), with 
z the transverse coordinate. In this expression, the subscripts label the phase of Xa 
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Figure 5: A schematic representation of the profile of the phase of the Xa fields within the 2- wall 
bound state in the upper plot. A'^ — 2 of the components have winding number 2/7V, while the 
remaining two have winding number 2/N — 1. In the lower plot the two 1-wall constituents are 
separated so that the fields all lie near the intermediate vacuum in the centre of the 2-wall. This 
configuration is unstable but can be used to extract the leading order interaction potential. 



in the relevant vacua, and d"' — x"" — x"'\ (o) are the deviations of each field from its 
value in the intermediate vacuum, while 7 is the phase of Zqi. These fields satisfy 
the linearised Bogomol'nyi equations given by 

N 

8,5" ^ e^'' Mf5\ where J]5" = 0, (4.11) 

a=l 

and = 5'"'^5c9(,VV(o) is the (complex) mass matrix in the intermediate vacuum, 
with Qab the Kahler metric. The additional constraint on the perturbations follows 
from f]^ = 1. 

The Bogomol'nyi equations decouple in the mass-eigenstate basis for the inter- 
mediate vacuum, where M is diagonal. In other words, since M is symmetric, the 
eigenvalues of the following system [40] 

are real and paired {m^, — m^}, for /j, ~ 1, . . . ,N — 1, and moreover {m^} can be 
identified with the mass eigenvalues in the intermediate vacuum. We deduce that 
the complex phase of the perturbations 5 arises purely from the integration constant, 
and can be fixed via comparison with (4.11). 
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With explicit solutions for the linear perturbations in hand, it is convenient to 
first use the Bogomol'nyi equation to write the potential in the form 

V{r) = -2Re [5(_io) (-^o) • ^^(oi) {zo + r)]+---, where = 9,5", (4.13) 

and the inner product is that associated with the Kahler metric. On substituting 
the appropriate solutions, we can read off the leading order potential, 

V{r, A) = -2 sin (^^^ J2 ^Me""""' + • • • , (4-14) 

where uj{\) is the relative phase between the central charges of the two constituent 
1-walls, as defined in (3.9). 

Note that the constraint apparent in the fundamental weight basis in (4.11) is 
reflected in the mass spectrum. For example, in the limit A = 0, we flnd 

M^{X = 0) oc =^ m^(A = 0) oc sin^ ^, (4.15) 

where is the Cartan matrix for affine SU(A^). The low-lying mass eigenvalues 
then appear in pairs, and depend nontrivially on N. 

It is now clear that the presence of the corresponding 2- walls depends purely on 
the value of a; (A), which is kinematic in nature as one should expect, since disconti- 
nuities in the BPS spectrum are allowed only when the corresponding bound states 
are marginally stable, namely when a; (A) = 0. Near these submanifolds the potential 
reduces to 

V{r, A) = -uj{X) m^e"™"'' + • • • , (4.16) 

which is linear in u and thus a discontinuity in the spectrum should apparently 
occur on crossing the curve. Moreover, the potential applies to all components of the 
multiplet (4.8), and so it would seem that the entire multiplet is removed from the 
spectrum when a;(A) < 0. However, before reaching this conclusion, we first need to 
consider the effect of quantum corrections. 

4.3 Quantum corrections 

Before drawing conclusions regarding discontinuities in the spectrum, it is important 
to understand the status of the leading order potential that we have derived. In 
particular, we are interested in submanifolds in parameter space on which the bound 
state is marginally stable, where a; (A) = 0, and the arguments above suggest that 
we can equivalently define the submanifold via the relation V{r, A) = 0. However, it 
turns out that the latter relation is not stable to quantum corrections, and more work 
is required. In actual fact quantum corrections, associated with the fermionic content 
of the model, are very important and become increasingly so near the CMS [54]. 
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Fortunately, one finds that the quahtative insight drawn from V{r, A) is nonetheless 
correct [54]. Indeed, this must be the case on the grounds that the condition a;(A) = 
describing the CMS is kinematic and protected by supersymmetry. 

To proceed, we will compactify the spatial worldvolume coordinates on an addi- 
tional circle of sufficiently small radius L. The low energy description of the world- 
volume theory, i.e. for scales well below 1/L, then reduces to quantum mechanics 
on a onc-dimensional moduli space coordinatiscd by r. In this regime, the leading 
order potential (4.16) is corrected by the fact that walls are 1/2-BPS, and thus the 
(dimensionally reduced) worldvolume theory must possess A/'=2 worldline supersym- 
metry. Such quantum mechanical systems arc well understood [56]. If wc denote by 
{Qi,Q2} the supercharges which act trivially on the putative BPS bound state, i.e. 

{Qif = {Q2f = {RL){T - \Z\) « {RL)T, (4.17) 

we expect to find a bound state if this system exhibits a unique vacuum; the multiplet 
structure is then reproduced on tensoring this state with the free centre-of-mass 
sector. 

The generic structure of J\f— 2 SQM was first described by Witten [56] , and we 
can realise the algebra as follows in terms of a quantum mechanical superpotential 
>VQM(r), 



Ql = Vr02 - >V^M(r)<7l] , (4.18) 

where — {RL){1/Ti + 1/T2)~^ is the effective reduced mass. The Hamiltonian is 
then given by 

T^QM = M - \Z\ = ^ [4 + (W'^uir)? + cTsWQuir)] , (4.19) 

where the second term is the classical potential, and the final term is a quantum 
correction oi 0{h). 

We have not explicitly allowed any nontrivial corrections to the kinetic term. 
In general one would expect that massive exchanges will introduce such terms, i.e. 
the coefficient of tt^ should take the schematic form g„. = 1 + 0{e~^'^). However, 
by a rescaling, always be absorbed into the potential, and this is implicitly 

the basis used here since we will in practice match to the full bosonic potential 

v{t) = 9"xw^y,r. 

Identifying the quantum mechanical superpotential Wquir) in general requires 
a detailed analysis of the model at hand, as in [54] . However, we can easily extract 
the general structure near the CMS by comparison with the leading order classical 
potential in (4.16) [57]. Having compactified the theory on T^, we can augment 
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(4.16) with the leading constant term corresponding to the binding energy of the 
two constituents. Near the CMS, the bosonic potential then has the generic form, 

V{r, A) = -Mruj{\f - uj{X) m^e""^"'' + • • • . (4.20) 

Comparing (4.19) with (4.20), we find 

WQM(r, A) = a;(A)(M,r) + J]exp(-m^r) + • • • , (4.21) 

and thus near the CMS, 

VQM(r, A) = ^Mru{\f - J] 

The additional quantum correction to the leading term becomes important near the 
CMS where it remains finite. 

Despite this correction, the classical equilibrium position survives as the maxi- 
mum of the ground state wavefunction, which is easily determined from (4.18), 

l^o) = exp(->VQM(r))|-) '-^ exp(-a;(A)M,r)|-) + • • • (4.23) 

where |±) are the eigenvectors of (73. In the limit r ^ 00 we have extracted only the 
leading exponential behaviour. It is apparent that this wavefunction is normalisable 
on only one side of the CMS. namely for > 0, which determines the existence 
domain for the BPS bound state. The intuition regarding discontinuities in the 
spectrum drawn from (4.16) therefore survives at the quantum level, although the 
details are somewhat different [54, 57]. 

4.4 Discontinuities in the BPS spectrum 

The analysis of the leading order potential allows us to reduce questions about the 2- 
wall multiplicity to purely kinematic issues concerning the dependence of the central 
charges on A. In particular, as discussed in section 3, discontinuities only occur on 
co-dimension one submanifolds of the parameter space - curves of marginal stability 
- where — 0. The analysis of this section has verified that such discontinuities 
do indeed occur in this case. 

Following our discussion in section 3. it is convenient for illustrative purposes to 
take large, so as to resolve the branch structure of the central charges as functions 
over the A-plane. To leading order in 1/iV, we recall that the angle uj takes the simple 
form 

-(A) = ^ 
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Figure 6: A plot of the circular (large N) CMS curve in the A-plane for the 2- wall composite. 
The 2-walls are absent in the shaded interior region. 

and one readily verifies that u is positive, and thus the BPS 2-wall bound states 
indeed exist outside the closed circle defined by |6A — 5p = 1, while the potential 
becomes repulsive, and the BPS bound states disappear from the spectrum in the 
interior shaded region of Fig. 6. 

In concluding this section, we will attempt to draw one further conclusion from 
these arguments relevant to the A ~ 1 regime, although this is strictly outside the 
regime of validity of the dynamical approach being used. Indeed, it is apparent that 
the expression for the relative phase at large N is singular at A = 1. This is the 
critical point at which the intermediate vacuum Wq collides with the corresponding 
confining vacuum Wn on the second branch, signalling the presence of additional 
light states. Moreover, for A^-even, A = 1 is also a Seiberg-Witten singularity and 
denotes a point where the confining branch intersects the Coulomb branch. 

For the composite 2-walls considered above, this suggests the possibility of form- 
ing a sandwich-like configuration where the outer regions are in confining phases - 
each associated in the near M = 2 regime with condensation of charge-one dyons 
and anti-dyons respectively - while the interior domain, between the two constituent 
1-walls, is in a massless phase at A = 1. Such configurations would be of interest 
as there are then generic arguments for the formation of fiux-tubes [34], i.e. open 
strings ending on the wall, and other features reminiscent of D-branes (see e.g. [58]). 
However, although not strictly valid for A = 1, the arguments above suggest that, at 
least within this system, there is no such BPS configuration, i.e. as we tune A ^ 1 
along a trajectory outside, but close to, the CMS the composite delocalises and on 
reaching A = 1 only the massless phase remains. 

5. Discussion 

We have presented a limited exploration of the spectrum of BPS domain walls within 
M = 1 gauge theories with adjoint matter, and specifically those characterised by 
a space of parameters associated with a classical polynomial superpotential. We 
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have been far from exhaustive and focused on a small subset of states which in an 
appropriate limit reduce to the BPS walls present in M= 1 SYM. The aim was simply 
to exhibit some of the features of the spectrum which, through the interpretation 
of the wall multiplicity in terms of the CFIV index, provides 'chiral data' on the 
parameter space. 

In this final section, we will return to some of the original motivations noted in 
Section 1, and comment on some extensions. 

5.1 Embedding within A/^=l* 

One motivation for studying the BPS spectrum was to gain insight into additional 
symmetries (or dualities) acting on the parameter space. With the aim of gaining 
a clearer picture of the quantum symmetries preserved by the BPS spectrum, it is 
helpful to consider the natural UV completion of the theory studied here within 
M= 4 SYM, a theory which exhibits S'-duality. Perturbing M= 4 SYM via the ad- 
dition of a cubic superpotential for the three adjoint chiral fields formally preserves 
modular invariance [59, 60] if we assign the mass terms modular weight (—5/6, 1/6), 
and the cubic couplings weight (—1,0). One may then construct the low energy 
effective superpotential by compactifying on x [11], and using the corre- 
spondence with the elliptic Calogero-Moser integrablc system [61, 10]. In partic- 
ular, the gauge invariant monomials are again identified with the action variables, 

(tr($P)) < > tr(LP), where L is the corresponding Lax matrix [62]. 

Suppressing the details, we can follow a similar procedure to that presented in 
Section 2 to determine the massive vacua. Restricting to those which correspond to 
a classically unbroken U(iV) gauge group, the central charges have the form [63], 

Z., = 2Y,nfm = ^Aki [(1 - - (1 - m^''] , (5.1) 

where 

e = 8/ — , (5.2) 

in terms of the mass parameters for the three chiral fields mj = {m, M, M), with 
m <^ M, and 

where {p, q) = (1, N), with A; = 0, . . . , 2N — 1 for the confining vacua, and (p, q) = 
{N, 1), with k — 0,N for the Higgs vacuum. E2 is the 2nd Eisenstein series, and 
C(r) is a vacuum-independent constant. With an appropriate choice of C(t), which 
corresponds to a particular basis for the action variables of the Calogero-Moser sys- 
tem, the combination ^X^ is modular invariant up to permutations of the vacua, and 
so the central charges have weight (—1/2, 1/2). Consequently, the wall tensions are 
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invariant up to permutation, and SL(2,Z) thus acts via permutations on the wall 
spectrum. 

The parameter space in this theory is now two-dimensional, coordinatised (for the 
massive U(A'") vacua) by {^, r}. One can show that in the decouphng hmit, M ^ oo, 
the spectrum of 2-walls interpolating between confining vacua exhibits the same CMS 
curve as found above, but whose profile is (not surprisingly) corrected by fractional 
instanton effects. The permutation symmetry of the vacua under rotations of A is now 
seen to follow directly from the symmetry under T-translations, namely r — > t + 1, in 
SL(2,Z) . A more complete analysis of the spectrum and symmetries, including the 
'classical' walls connecting confining and Higgs vacua, will appear separately [63]. 

5.2 On the anomaly multiplet 

Another motivation for this work was to try and understand whether the matrix 
model approach could be used to extract the chiral data associated with the BPS 
spectrum - namely to compute the CFIV index - in a purely four-dimensional con- 
text. An important point to recognise is that, while the confining vacua are described 
by 1-cut solutions to the matrix model, this is insufficient to describe BPS states. 
Indeed, the field profiles required to construct domain walls will necessarily pass 
through generic regions in moduli space where, in the A/'=2 limit, the gauge group is 
maximally broken. This implies that we require knowledge of the large M solution 
to the matrix model with the maximal number of cuts. On reflection, this should 
not be too surprising on comparison with the approach used here, where we made 
use of the full N"— 2 Seiberg-Witten curve. 

One may then enquire as to how the nontrivial wall multiplicity, described in 
section 4 [53, 49], arises in the pure A/'= 1 SYM limit where, having integrated 
out the adjoint field, the effective supcrpotential W = VV(S') depends only on a 
single gluino condensate field S. Even discounting non-analyticities, a supcrpotential 
depending on a single field can at most describe unique interpolating solutions, and 
this is insufficient to explain the multiplicities found in section 4. One may then 
ask how the information about this multiplicity is encoded in four dimensions, if 
at all? To this end, it is worth recalling the resolution of a (perhaps) analogous 
puzzle that arises on compactification to 1+lD on a 2-torus (see also [17]). Since 
we are computing the index trF(— 1)^, one might anticipate that such an additional 
reduction should not affect the conclusions. 

We can ignore the presence of the adjoint chiral field, so the reduced supcrpo- 
tential follows directly from a reduction of the 3D affine-Toda supcrpotential. 



where Xa — e with A2D = 47r^(it!L)A^^]^, and we have inserted a Lagrange 
multipher to exphcitly enforce the decoupling of the overall U(l) factor. One can 




(5.4) 
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recognise this as the 1+lD mirror [64] of the N'= 2 CP^~^ sigma model where we 
interpret S as the U(l) field strength within the corresponding linear sigma model. 
This relationship is also consistent with a direct reduction of pure SYM, since the 
moduli space of flat SU(iV) connections on a torus is also CP^~^ [65]. 

If we now integrate out E from (5.4), we simply enforce the decoupling of the 
central U(l), and we can reproduce the vacuum structure and nontrivial kink mul- 
tiplicity in the same manner as described in section 4 [50]. However, we could also 
choose to integrate out the v^s, or equivalently, the homogeneous CP^~^ chiral 
fields. From (5.4), we can do this at tree level to find. 



This reduced theory describes, as it must, the same vacuum structure but, due to 
Gauss law, actually cannot describe any nontrivial kinks as one must satisfy the 
constraint S|_oo = '^\+oo [66]. As noted by Witten, the resolution of this puzzle 
is that the kinks are actually charged under the additional CP^~^ fields that have 
been integrated out, and one must account for the coupling of E to the corresponding 
current [66, 67]. In this sense, the superpotential (5.5), while sufficient for describing 
the vacua, is insufficient to fully describe the BPS spectrum. 

Although there is no direct analogue of the Gauss law constraint in 3+ ID, where 
we lift (5.5) to the Veneziano-Yankiclowicz superpotential [38], one could make quite 
a close analogy if we were to interpret 5* as the field strength of a linear multiplet 
[68] containing a 3-form. While we will not pursue this further here, this issue is 
intriguing as these kinks are precisely the dimensional reduction of the BPS wall 
configurations that we have been counting. 
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